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7n wiro velocities of several cell types have been measured using computer controlled video mi- 
croscopy, which allowed to record the cells' trajectories over several days. On the basis of our large 
data sets we show that the locomotion activity displays a universal exponential distribution. Thus, 
motion resulting from complex cellular processes can be well described by an unexpected, but very 
simple distribution function. A simple phenomenological model based on the interaction of various 
cellular processes and finite ATP production rate is proposed to explain these experimental results. 



Tissue cell migration |TJ has been studied in a num- 
ber of biological processes, including the development of 
the nervous system H, immune reactions [0 and tumor 
spreading m. Beside the inherent biological importance, 
cell locomotion plays a key role in a number of interesting 
large scale pattern formation phenomena like cell sorting 
p-M, or aggregation of amoebae H. Although the be- 
havior of swimming bacteria has been well studied H and 
served as a basis to understand collective microbiological 
phenomena |1C| |, rather little is known about the basic 
statistical features of eukaryotic cell migration on time 
scales comparable to the time of the pattern formation 
processes. 

The locomotion of cells emerges from parallel cyclic 
subprocesses including cell-attachment to and detach- 
ment off the substrate, growth and collapse of filopo- 
dia/lamellopodia, the dynamic maintenance of appropri- 
ate cell surface composition and displacement of the cell 
mass. This complexity suggests two consequences. On 
one hand, strong correlation cannot be expected between 
the actual cell velocity and the concentration of a single 
(or a few) cell biological regulatory factor(s) fll|j . On 
the other hand, interacting complex units often results a 
characteristic macroscopic behavior which can be under- 
stood by mathematical models [p|-p|p0| , thus some sort of 
statistical regularities might be expected in the cellular 
behavior as well. 

In order to obtain data providing a basis for an exten- 
sive statistical analysis we have developed a computer 
controlled system which automatically collects data and 
maintains the cell culture conditions in a chamber at- 
tached to the microscope for several days. We found 
that the locomotion activity displays an exponential dis- 
tribution, i.e., in spite of the complex cellular processes, 
the motion of a cell can be well described by a simple 
universal distribution function. This observation can be 
explained in terms of a phenomenological model based 
on the competition of various cellular processes for the 
finite free energy (ATP) resources. 

With the aim of applying statistical approaches in the 
study of cellular locomotory activity, long-term cell mi- 



gration patterns were recorded in low-density monolayer 
cultures (see Table 1). The cellular trajectories (Fig. 1) 
were tracked using a computer controlled phase contrast 
videomicroscope system. Due to a high degree of autom- 
atization we could collect and analyze a large set of data 
(about 100,000 cell positions), significantly exceeding the 
amount of information evaluated in previous works on 
ceU motility IggHHH. 

Average velocities were calculated over a time period 
of one hour (t ) asuj(i) = |.f i (i + t /2) — Xi{t — t /2)\/t , 
where the location of the ith cell at time t is denoted by 
Xi{t). Remarkable fluctuations were found in the loco- 
motion activity of individual cells, a phenomenon which 
has also been reported in living slices of developing cortex 
[g| , in vitro migration of neurons |2(| and fibroblasts [|15| . 
The characteristic time scale of the fluctuations was esti- 
mated to be a few hours, a value similar to that reported 
for fibroblasts JlEJ ]. 

Due to the low cell density, most cells migrated freely, 
but the effects of various cell interactions could be also 
investigated. Since it is known that high cell density can 
significantly decrease the cell locomotion activity pi[ , 
we investigated whether the observed velocity fluctua- 
tions were caused by changes in local intercellular connec- 
tions. Cell to cell contacts were identified on the images 
and the locomotion velocities Vi(t) were sorted into two 
sets ("solitary" / "interacting") depending on whether 
the cell in question had or had not visually observable 
contacts with other cells. Cumulative velocity distribu- 
tion functions F(v) were calculated for time periods of 
20 hours giving the probability of the relation Vi (t) < v 
to hold for a randomly selected cell i and a random time 
t. Comparing F(v) of "solitary" cells to that of "inter- 
acting" cells (Fig. 2a), two almost overlapping ranges 
of velocities (0 — 80 /xm/h and — 100 /im/h, respec- 
tively) were found. It is also known |j[] that during migra- 
tion cells leave various extracellular molecules and other 
polipeptides attached to the substrates. To elucidate the 
effect of these cell trails on the velocity fluctuations we fil- 
tered our database in such a way that we kept only those 
cell positions which distance was greater than 20^m (i.e., 



twice the typical cell size) of the trajectory of other cells. 
Again, no qualitative difference was found in the velocity 
fluctuations suggesting that they are inherent property of 
the cells. 

One of our main findings is that F(v) of all cultures 
investigated could be very well fitted by an exponential 
distribution 
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within measurement errors in the entire velocity range 
studied (Fig. 2b). This result indicates that the ex- 
ponential velocity distribution is likely to be a general 
characteristic feature of in vitro cell motility for a broad 
class of surface attached cells. 

The above velocity distributions were found for rela- 
tively large (ss 100 cells) populations. On the level of 
individual cells the exponential behaviour (|lj) can be in- 
terpreted in two ways: (i) The culture is inhomogeneous, 
i.e., slower and faster cells can be distinguished on the 
bases of well preserved phenotypic properties. In this 
case the exponential F(v) distribution can reflect the ra- 
tio of the slow and fast cells in the culture, while the 
velocity fluctuations of the individual cells can show an 
arbitrary distribution, (ii) If the culture is homogeneous, 
then almost all cells exhibit the same distribution of ve- 
locity fluctuations, i.e., Fi(v) ~ F(v) holds for each cell i, 
where Fi (v) denotes the distribution function of the cell 
i. In this scenario the average velocity (vi(t)) of each cell 
would be the same if we could calculate the time averages 
over an infinitely long time. Since the time averages are 
calculated over a finite time r only, for the distribution 
of the average velocities we can expect a Gamma distri- 
bution with a parameter s = r/to, where the correlation 
time of the process is denoted by t [g3| . 

To decide between the above alternatives Fi(v) and 
the average velocity (vi(t)) over a period of r = 16h were 
calculated for each cell i. Fig. 3a and 3b show that the 
experimental data clearly support the homogeneous al- 
ternative. As an example, the data shown in Fig. 3b 
can be well fitted by a normalized Gamma distribution 
Fr(u,s) = L dzs s z s ~ 1 e~ sz /(s — 1)! with a parameter 
value of s = 9. The fitted value of s is consistent with 
the finding that the correlation time scale of the velocity 
fluctuations is in the order of hours. 

The fact that the distribution of cellular velocities fol- 
lows a simple exponential function is an unexpected find- 
ing. Several subprocesses of cell locomotion (receptor 
binding, membrane exocitosis, lamelopodium/filopodium 
formation, etc.) display either Poisson distribution |y| 
or periodic behavior |2^] if studied on time scales short 
enough to exclude non-stationarities of the cell state. 
However, the fluctuations of cell velocities recorded in our 
experiments were observed on significantly longer time 
scales than those listed above. Since in a few hours many 
lamellopodia are formed and many vesicles are fused to 
the membrane, etc., the relative frequency fluctuation 
of such stationary processes must be small and close to 



Gaussian. In agreement with this expectation, previous 
works on cell locomotion - although not calculating F(v) 
explicitly - have set up phenomenological models pre- 
dicting Gaussian velocity distribution |l^Jlq , [l8|] . Also, 
a close to Gaussian distribution was found in |6j where 
motion within a cell aggregate was studied. In that case, 
however, there was a strong interaction between the cells 
and the displacement of a single cell was also a result of 
the activity of the neighboring cells. 

Within the time scale of the velocity fluctuations both 
the cellular environment and the pattern of gene expres- 
sion in the cells can be considered as stationary. This 
time, on the other hand, is long enough to permit changes 
in cell motility as a consequence of changes in the concen- 
tration of several intracellular regulatory factors. In the 
following we demonstrate with a very simple cell model 
that the exponential distribution of cell velocity fluc- 
tuations can indeed reflect such stochastic intracellular 
changes. 

Since the cell locomotion processes are cyclic and dis- 
sipative, the observed velocity is proportional to the 
amount Jq of free energy dissipated (number of ATP 
molecules hydrolized) in the appropriate chemical reac- 
tions contributing to cell motility: 



v(t) = M)(i), 



(2) 



where fi denotes the conversion factor between velocity 
and energy dissipation. It is natural to assume that fj, is 
constant, i.e., the subprocesses are synchronized in such a 
manner that a close to optimal distribution of the energy 
inflow is maintained. 

The individual reaction rates are determined on one 
hand by the ATP concentration c and on the other hand 
by the concentrations of various regulatory factors. Some 
of these factors maintain the motility efficiency fj,, while 
other factors are products of different cellular processes, 
with the potential of altering the cell locomotion activ- 
ity. These latter factors can change independently from 
each other: while some of the changes increase, some oth- 
ers decrease the cell motility. Thus, we assume that the 
temporal changes of Jq can be represented as 
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dt 
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where go(t) > and £,o(t) is a stochastic variable with 
zero mean and Gaussian distribution. 

Finally, we assume that the cellular ATP production 
rate P is limited by the ATP production capacity of the 
cell, so the various cellular processes (among them cell 
locomotion, protein synthesis, DNA duplication, etc.) 
compete for the ATP available. Since the changes of c 
can be calculated as a difference in the production and 
the consumption rates, the equations 
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0,1,...,7V (4) 



fc=0 



are obtained, where Ji denotes the rate of free energy 
dissipation in a given cellular process C, while gi and ^ 
are analog quantities to go and £o> respectively. 

It can be seen that Eqs. (|J) describe a diffusion pro- 
cess in the vicinity of the J2k=o Jk = P hyperplane. Ac- 
cordingly, the probability distribution of the "cell state" 
J = (Jo, J\, ..., Jjv) on this hyperplane is uniform. Inte- 

fiting the distribution of J and taking into account Eq. 
we obtain the experimentally observed exponential 
behaviour 
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where p = P/(N + 1) is the average free energy consump- 
tion rate of a given cellular process. 

There is a formal similarity between the distribution 
given by Eq. (g) and the Boltzmann energy distribution 
within systems in thermal equilibrium. In the latter case, 
the uniform micro canonical distribution can be derived 
using either quantum mechanics or Liouville's theorem 
for Hamiltonian systems while here it is an interesting 
feature of the interacting regulatory processes. 

Although Eq. (Q) is certainly oversimplified and the 
free energy dissipation rates Ji are not identified with 
actual biochemical reactions, the exponential behavior 
of Eq. (g|) is expected to hold for a wide class of models 
(e.g., non-constant P, considering explicitly the interac- 
tions of regulatory networks, etc.). Our simple calcula- 
tion also predicts that exponentially distributed fluctua- 
tions can be expected in numerous other cellular activi- 
ties as well. 

We thank Z. Csahok and O. Haiman for discussions, 
Zs. Kornyei and B. Hegediis for help in the data analysis, 
I. Fazekas for providing us the human glioblastoma cell 
lines. Supporting grants OTKA T019299, F017246 and 
FKFP 0203/1997 were substantial for our work. 



[2] 



[9 
[10 



[11 
[12 

[is; 

[14 
[15. 

[is; 

[17 
[18 
[19 

[20 
[21 

[22 



[23 



[1] For review see, e.g., Dynamics of Cell and Tis- 
sue Motion, ed. W. Alt, A. Deutsch and G. Dunn 



(Birkhauser, Basel, 1997); D.A. Lauffenburger, and A.F. 

Horwitz, Cell 84, 359 (1996). 

N.A. O'Rourke et al. , Science 258, 299 (1992); J. Fok- 

Seang et al. , Dev. Biol. 171, 1 (1995); C.E. Krull et al. , 

Development 121, 3733 (1995). 

P. A. Negulescu et al. , Immunity 4, 421 (1996). 

M.R Chicoine and D.L. Silbergeld, J. Neurosurg. 82, 615 

(1995). 

R.A. Foty et al. , Phys. Rev. Lett. 72, 2298 (1994). 

J.C.M. Mombach and JA. Glazier, Phys. Rev. Lett. 76, 

3032 (1996). 

JP. Rieu et al. , Phys. Rev. E. 57, 924 (1998). 

D.A. Kessler and H. Levine, Phys. Rev. E. 48 4801 

(1993); B.N. Vasiev et al. Phys. Rev. Lett. 73 3173 

(1994); Hofer et al. , Physica D 85 425 (1995). 

H.C. Berg and E.M. Purcell, Biophys. J. 20, 193 (1977). 

E.O. Budrene and H.C. Berg, Nature 376, 49 (1995); 

J.O. Kessler, Nature 313, 218 (1985); M. Matsushita and 

H. Fujikawa, Physica A 168, 498 (1990); E. Ben- Jacob 

et al. Nature 386, 46 (1994). 

J.T.H. Mandeville et al. , Biophys. J. 68, 1207 (1995). 

K. Schlett and E. Madarasz, J. Neurosci. Res. 47, 405 

(1997). 

W. Twal et al. , in preparation 

E. Madarasz et al. , Neuroscience 43, 211 (1991). 

G.A. Dunn and A.F. Brown, J. Cell. Sci. Suppl. 8, 81 

(1987). 

M.R. Parkhurst and W.M. Saltzman, Biophys. J. 61, 306 

(1992). 

G. Fishell et al. , Nature 362, 636 (1993). 

H.M. Buettner et al. , Dev. Biol. 163, 407 (1994). 

H. Hu et al. , Neuron 16, 735 (1996). 

H. Wichterle et al. , Neuron 18, 779 (1997). 

see, e.g., The molecular biology of the cell, second edition, 

ed. B. Alberts et al. (Garland, New York, 1989) pp 673. 

B. Hintz and O. Brostenau, in Dynamics of Cell and 

Tissue Motion, ed. W. Alt, A. Deutsch and G. Dunn 

(Birkhauser, Basel, 1997), pp 15. 

The time averaged velocity can be written as (vi(t)) T = 

ELo*(*)A = (V*)E: = i(lAo)Ei:ttT°«(*). where 
s = r/£o- Since the terms in the second sum are posi- 
tively correlated, we estimate the average by (vi(t)) r ~ 
(V s )El=i"( fl ^)' which expression is a sum of s inde- 
pendent, exponentially distributed random variable. 



FIG. 1. Cellular trajectories of NE-4C cells in a field of 1110 /imx 840 /xm, during the first 20 hours in culture. Microscopic 
images were taken with 10 x objective magnification every 3-4 minutes. The cells in each 3rd image were tracked manually 
with a precision of ca. 5/xm, which is comparable with the size of the cells (10 (im). 



FIG. 2. (a) Cumulative velocity distribution functions F(v) calculated from the trajectories of Fig. 1. for three subsets of 
cells: "interacting", "solitary" and cells which do not cross the trails left on the substrate (see text for details). The similarity 
of the curves indicates that the fluctuations are an inherent property of the cells. The horizontal bars correspond to the error 
of velocity determination due mainly to the method of cell positioning. The vertical bars mark the systematic error ( w 5% ) 
of F(v), which estimation was based on the difference between the distribution functions calculated for the first and the second 
10 hours of the record. The solid line is a fitted exponential distribution function. The inset shows the same distribution 
functions on a linear- logarithmic plot, (b) The velocities u = v/{vi(t)) normalized by the average velocity of the population 
show an exponential distribution in all the 13 cultures investigated, demonstrated by the linear region covering four decades in 
the linear-logarithmic plot. The figure shows a typical result for each cell line. 

FIG. 3. (a) The average function of the single-cell velocity distributions Fi(v) = y^ n _ Fi(v)/n show the same exponential 
behaviour as the cumulative distribution function F(v) of the whole cell population shown in Fig.l. The dotted lines are 
an exponential fit of F(v) in the corresponding velocity regime, (b) The distribution function of the normalized average 
cellular velocities (ui) T — J Vi(t)dt/(r{ Vi(t)}) is presented, both as linear and lin-log plots. The solid line is a fitted Gamma 
distribution. 



TABLE I. The most important features of the cultures investigated 



Culture 


reference 


cell density a 
[cells/mm 2 ] 


duration 
[h] 


positions 


(Vi(t)) 

[/im/h] 


mouse neuroepithelial progenitor (NE-4C) 


[12] 


20, 40, 200 


24, 24, 12.5 


2000, 4000, 12500 


6.7, 8.7, 17 


mouse fibroblast (NIH-3T3) 


ATCC' 


170 


20 


5000 


12.1 


human glioblastoma (HA) 


OITP 


48, 65 


50, 50 


13500, 19500 


11.1, 11.6 


human glioblastoma (HC) 


OITI 


126, 80 


50, 50 


37800, 24000 


6.5, 7.0 


human breast adenocarcinoma (MDA MB231) 


ATCC, [13] 


40 


24 


2500 


5 


rat primary astroglia 


[14] 


200 d 


50, 50, 36, 30 


3500, 3500, 3000, 2500 


8 



a at seeding; b American Type Culture Collection; c National Institute of Neurosurgery, Hungary; ^migration from cell aggregates 
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